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Static properties of multiskyrmions with baryon numbers up to 8 are calculated, including 
momenta of inertia and sigma-term. The calculations are based on the recently suggested 
SU (2) rational map ansaetze. The spectra of baryonic systems with strangeness, charm 
Q\ \ and bottom are considered within a "rigid oscillator" version of the bound state soliton 

OS i model. The binding energies estimates are made for the states with largest isospin which 

can appear as negatively charged nuclear fragments, as well as for states with zero isospin 
- light fragments of "flavoured" nuclear matter. Our results confirm the previously made 
observation that baryonic systems with charm or bottom quantum numbers have more 
. chance to be stable with respect to strong interactions than strange baryonic systems. 
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1 Introduction 



The topological soliton models, and the Skyrme model among them [fy, are attractive 
because of their simplicity and the possibility that they may describe well various prop- 
erties of low energy baryons. The models of this kind provide also a very good framework 
within which to investigate the possibility of the existence of nuclear matter fragments 
with unusual properties, such as flavour being different from u and d quarks. In addition to 
OS ' being important by itself, this issue can have important consequences in astrophysics and 
cosmology. In particular, the formation and subsequent decay of such fragments could be 
Ph! important in the early stages of the evolution of the Universe. It is well known that the 
relativistic many-body problems cannot be solved directly using the existing methods, 
and the chiral soliton approach may allow us to overcome some of these difficulties. 

The description of skyrmions with large baryon numbers has been perceived as 
being complicated because the explicit form of the fields has been not known. A recent 
observation that the fields of the SU(2) skyrmions can be approximated accurately 
by rational map ansaetze giving the values of masses close to their precise values, has 
simplified considerably the task of such studies. Similar ansaetze have also been recently 
presented for SU(N) skyrmions (which are not embeddings of SU(2) fields) [[5]. 

In this paper we use the SU(2) rational map ansaetze as the starting points for the 
calculation of static properties of bound states of skyrmions necessary for their quanti- 
zation in the SU(3) collective coordinate space. The energy density of the 5 = 3 config- 
uration has tetrahedral symmetry, of B = 4 - the octahedral (cubic) one [f|], of B = 5 - 
-D 2 d-symmetry, of B = 6 - £> 4d , of B = 7 - dodecahedral symmetry, and of B = 8 - D 6d - 
symmetry || §], etc. 

The minimization, with the help of a 3-dimensional variational program || , lowers 
the energies of these configurations by few hundreds of Mev and shows that they become 
local minima in the SU(3) configuration space. The knowledge of the so-called "flavour" 
moment of inertia and the E-term allows us then to estimate the flavour excitation ener- 
gies. The mass splittings of the lowest states with different values of strangeness, charm 
or bottom are calculated within the rigid oscillator version of the bound state approach. 
The binding energies of baryonic systems (BS) with different values of flavours are also 
estimated. 

To reduce theoretical uncertainties we consider the differencies between the binding 



energies of BS with flavour F and the ground state for each value of B. These ground 
states are the deuteron for B = 2, the isodoublet 3 H- 3 He for B = 3, 4 ife for £? = 4, etc. 
These differencies, being free of many uncertainties, in particular of the poorly known loop 
corrections to classical masses, show the tendency of the flavoured BS to be more bound 
than the (u,d) ground states (for heavy flavours), or to be less bound, as for strangeness. 
Of course, it is an assumption that the ground states of multiskyrmions correspond to 
ordinary nuclei. However, this is a natural assumption if we believe that effective field 
theories describe nature. 

In the next Section the static characteristics of multiskyrmions are described. 
Flavour excitation energies and zero mode corrections to the energies of multibaryons are 
considered in Section 3. Section 4 contains estimates for the binding energies of baryonic 
systems with different values of flavours, and our conclusions are given in Section 5. 



2 Static characteristics of multiskyrmions 

We consider here simple SU(3) extensions of the Skyrme model 0: we start with SU(2) 
skyrmions (with flavour corresponding to (u,d) quarks) and extend them to various SU(S) 
groups, such as, (u,d,s), (u,d,c), or (u,d,b). 

We take the Lagrangian of the Skyrme model which, in its well known form, 
depends on parameters F v , F D and e and can be written in the following way 0: 

C = ^TrlJ" + ^Tr[l M , l v f + ^^Tr(U + - 2)+ 

+ F v m D- F * m * Tr(l - V3\ 8 )(U + f/t _ 2 ) + Fl3 ~ F * Tr(l - V3X 8 )(Ul^ + l^tf), (1) 
24 4o 

Here U e SU(3) is a unitary matrix incorporating chiral (meson) fields, and = U^d^U. In 
this model F^ is fixed at the physical value: F^ = 186 Mev and M D is the mass of K, D or 
B meson. The ratios Fd/F^ are known to be 1.22 and 1.7 ±0.2 for, respectively, kaons and 
D-mesons. 

The flavour symmetry breaking (FSB) in the Lagrangian is of the usual form, and 
is sufficient to describe the mass splittings of the octet and decuplets of baryons within 
the collective coordinate quantization approach @. 

The Wess-Zumino term, to be added to the action, which can be written as a 
5-dimensional differential form || plays an important role in the quantization procedure: 

SWZ = 24^ L d^e^Tr^lM^), (2) 

where il is a 5-dimensional region with the 4-dimensional space-time as its boundary and 
= U^dfjJJ. Action (2) defines important topological properties of skyrmions, but it does 
not contribute to the static masses of classical configurations || 0. Variation of this 
action can be presented as a well defined contribution to the Lagrangian (integral over 
the 4-dimensional space-time). 

We begin our calculations, however, with U e SU(2). The classical mass of SU(2) 
solitons, in the most general case, depends on 3 profile functions: /, a and /3 and is given 
by 

M cl = | j £ji $ + £ + j*] + _L [lijtf + { l 2 T 3 f + [ytf] + -Flml(\ - C/ )|d 3 r (3) 

Here l k are the SU(2) chiral derivatives defined by WdU = UkTk, where k = 1, 2, 3. The 
general parametrization of U for an SU(2) soliton we use here is given by Uq = c/ + s/fn 



with n, = 



s a s/3, sf = sinf, cf = cosf, etc. For the rational map ansatz we 



are using here as starting configurations, 

2ReR(£) _ 2ImR{£) 



i-TOI s 



where R(£) is a ratio of polynomials of the maximal power B of the variable £ = tg(0/2)exp(i(j>), 
9 and <p being polar and azimuthal angles defining the direction of the radius-vector f. 
The explicit form of R(£) is given in for different values of B. 

The "flavour" moment of inertia plays a very important role in the procedure of 
1, see formulas (9), (10) below, and for arbitrary SU{2) skyrmions 



SU(3) quantization 
is given by [[H], [T!| : 



Qf = IJ(1- c f ) [F* + i((9/) 2 + s%daf + S y a 0pf)]d 3 f. 



It is simply connected with Qy of the flavour symmetric case (F D 

e F = ef + {fI/f* i)r/4, 



(4a) 



(46) 



with r defined in (5) below. The isotopic momenta of inertia are the components of the 
corresponding tensor of inertia. They have been discussed in many papers, see e.g. [01- 
||12|| , so, we will not present them here. For majority of multiskyrmions we discuss, this 
tensor of inertia is close to the unit matrix multiplied by the isotopic moment of inertia ©t- 
This is exactly the case for 5 = 1 and, to within a good accuracy, for 5 = 3,7. Considerable 
deviations take place for the torus with 5 = 2, and smaller ones for B = 4, 5, 6 and 8, see 
Table 1. The quantity r (or E-term), which defines the contribution of the mass term 
to the classical mass of solitons, and f are used directly in the quantization procedure. 
They are given by: 
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(l- C/ )d 3 f, 



Cf[(df? 



\(da) 



The following relation can also be established: f = 2(M^/F 2 



4 



(5) 



where is the 



second-order term contribution to the classical mass of the soliton, and 0f, fe is the Skyrme 
term contribution to the flavour moment of inertia. The calculated masses of solitons, 
momenta of inertia Qf, ©t, r or E-term and f are presented in Table 1. 



B 


Mci 


e (0) 


©T 


9 T 3 


r 


f 




UJ C 




Cs 


c c 


Cb 


C s 


c c 


1 


1.702 


2.04 


5.55 


5.55 


4.83 


15.6 


0.309 


1.542 


4.82 


0.28 


0.27 


0.52 


0.54 


0.91 


2 


3.26 


4.18 


11.5 


7.38 


9.35 


22 


0.293 


1.511 


4.76 


0.27 


0.24 


0.49 


0.53 


0.90 


3 


4.80 


6.34 


14.4 


14.4 


14.0 


27 


0.289 


1.504 


4.75 


0.40 


0.37 


0.58 


0.60 


0.92 


4 


6.20 


8.27 


16.8 


20.3 


18.0 


31 


0.283 


1.493 


4.74 


0.47 


0.44 


0.62 


0.64 


0.92 


5 


7.78 


10.8 


23.5 


19.5 


23.8 


35 


0.287 


1.505 


4.75 


0.42 


0.40 


0.60 


0.62 


0.92 


6 


9.24 


13.1 


25.4 


27.7 


29.0 


38 


0.287 


1.504 


4.75 


0.48 


0.46 


0.63 


0.67 


0.93 


7 


10.6 


14.7 


28.9 


28.9 


32.3 


44 


0.282 


1.497 


4.75 


0.48 


0.46 


0.64 


0.66 


0.93 


8 


12.2 


17.4 


33.4 


31.4 


38.9 


47 


0.288 


1.510 


4.79 


0.49 


0.47 


0.64 


0.67 


0.93 



Table 1. Characteristics of the bound states of skyrmions with baryon numbers up to B = 8. The classical 
mass of solitons M c i is in Gev, momenta of inertia Of, ©t and ©t.3, T and T - in Gew -1 , the excitation 
frequencies for flavour F, u>f in GeV . c StC .b and c SiC are the hyperfine splitting constants for multibaryons 
defined in Eq.(21). The constant Cb is close to 0.99 for all B and is not included into the Table. The 
parameters of the model F„ = 186 MeV, e = 4.12. The accuracy of calculations is better than 1% for 
the masses and few % for other static characteristics of solitons. The B — 1 quantities as well as B = 2 



quantities for the torus calculated previously, are shown for comparison. 



As can be seen from Tables 1 and 2, there are two "islands" of stability for baryon 
numbers considered here: at B = 4, which is not unexpected, and for B = 7, and this is 
something new and unexpected. So far, this property seems to be specific to the Skyrme 
model. The difference between Q T and 9t,3 is maximal for the toroidal B = 2 configuration 
and decreases with increasing B. It vanishes for B = 3 and 7 precisely. The accuracy of 
calculations decreases with increasing baryon number. It is difficult to estimate it for such 
quantities as u f ,b and c FiB where it depends also on the particular method of calculation 
- the rigid oscillator model in our case. 

The behaviour of static properties of multiskyrmions and flavour excitation fre- 



quencies shown in Table 1 is similar to that obtained in [|22[ for toroidal configurations 
with B = 2,3,4. We note that the flavour inertia Qf.b increases with B almost proportion- 
ally to B. 



3 Flavour excitation frequencies and ~ l/N c zero mode correc- 
tions 

To quantize the solitons in their 5(7(3) configuration space, in the spirit of the bound state 
approach to the description of strangeness proposed in [|13|-[]14| and used in |l5fl-[|17lj, we 
consider the collective coordinate motion of the meson fields incorporated into the matrix 
U: 

U(r, t) = R(t)U (O(t)r)rf(t), R(t) = A(t)S(t), (6) 

where U is the SU(2) soliton embedded into SU(3) in the usual way (into the upper left 
hand corner), Ait) e SU(2) describes SU(2) rotations and S(t) e SU(3) describes rotations 
in the "strange", "charm" or "bottom" directions and 0(t) describes rigid rotations in 
real space. For definiteness we consider the extension of the (u, d) SU(2) Skyrme model 
in the (w, d, s) direction, when D is the field of if-mesons but it is clear that quite similar 
extensions can also be made in the directions of charm or bottom. So 

Sit) = exp(W(t)), V(t) = D a(t)X a , (7) 

a=4,...7 

where A a are Gell-Mann matrices of the iu,d,s), (u,d,c) or (u,d,b) SU(3) groups. The 
(u,d,c) and (u,d,b) SU{3) groups are quite analogous to the (u,d,s) one. For the (u,d,c) 
group a simple redefiniton of hypercharge should be made. For the (u,d,s) group, L> 4 = 
(K+ + K~)/V2, D 5 = i(K+ - K~)/V2, etc. For the (u,d,c) group D t = {D° + D°)/V2, etc. 

The angular velocities of the isospin rotations u are defined in the standard way 
0: At A = -iQf/2. We shall not consider here the usual space rotations in detail because 
the corresponding momenta of inertia for BS are much greater than the isospin momenta 
of inertia, and for the lowest possible values of angular momentum J, the corresponding 
quantum correction is either exactly zero (for even B), or small, see also formulas (17) 
and (21) below. 

The field D is small in magnitude. In fact, it is, at least, of order l/^fW c ^ where N c 
is the number of colours in QCD, see Eq. (14). Therefore, the expansion of the matrix S 
in D can be made safely. 

The mass term of the Lagrangian (1) can be calculated exactly, without expansion 
in the powers of the field D, because the matrix S is given by S = l—iV sind/d-V 2 (l-cosd)/d 2 



with d 2 = TrV 2 . We find that 



AC 



M 



FWd 



(8) 



The expansion of this term can be done easily up to any order in d. The comparison of 
this expression with AL M , within the collective coordinate approach of the quantization of 
SU{2) solitons in the SU(3) configuration space [[H|-[|l2|, allows us to establish the relation 
sin 2 d = sin 2 v, where v is the angle of the A4 rotation, or the rotation into the "strange" 
("charm", "bottom" ) direction. 

After some calculations we find that the Lagrangian of the model, to the lowest 
order in the field D, can be written as 



-M, 



cl,B 



F 2 

I r D 2 2 
r B ( ~j^ m D - m 



f B (F 2 D -F 2 ) 



D^D-i^-(D^b-t)^D). (9) 



Here and below D is the doublet K+, K° (D°, D~, or B+, B°): d 2 = TrV 2 = 2D^D. We have 
kept the standard notation for the moment of inertia of the rotation into the "flavour" 

0- 



12 1 ; different notations are used in |15|, [T(| (the index 



direction <d F for C , 9 b or 9 S 

denotes the charm quantum number, except in N c ) . The contribution proportional to f B 
is suppressed in comparison with the term ~ r by a small factor — (Fg - F 2 )/m 2 D , and is 
more important for strangeness. 

The term proportional to N C B in (9) arises from the Wess-Zumino term in the 
action and is responsible for the difference of the excitation energies of strangeness and 
antistrangeness (flavour and antiflavour in the general case) ||13||-fT6f . 



Following the canonical quantization procedure the Hamiltonian of the system, 
including the terms of the order of N°, takes the form [15|, 



M, 



cLB 



46 



-n f n 



F,B 



Bm D 



?b{F 2 d 



Ft 



N 2 B 



2 1 



169 



F,B 



D^D 



+ i 



N r B 



89 



F,B 



(D^U-UW), (10) 



where m 2 D = (F^/F 2 )m 2 D 



The momentum n is canonically conjugate to variable D 



(see Eq.(18) below). Eq. (10) describes an oscillator-type motion of the field D in the 
background formed by the (u, d) SU(2) soliton. After the diagonalization, which can be 
done explicitly following |15|, [T(J , the normal-ordered Hamiltonian can be written as 

(11) 



H, 



M cLB + uj F . B a)a + uJF.E^b + 0(l/N c ) 



with a^, fei being the operators of creation of strangeness (i.e., antikaons) and antis- 
trangeness (flavour and antiflavour) quantum number, uj f , b and luf,b being the frequencies 
of flavour (antiflavour) excitations. D and n are connected with a and b in the following 
way p], |TB[: 

1 



D l = 



y/N c Bn FiB 



IT 



y/N c BfJ, FiB 



2i 



with 



Hf,b = [1 + 16(fh 2 D T B + (F 2 - F 2 )f B )Q F , B /{N c B) 2 ] 1 l 2 . 
For the lowest states the values of D are small: 

-1/4 



D - WT B Q FB ffih+N^B 



(12) 



(13) 



(14) 



and increase, with increasing flavour number like (2|F| + l) 1 / 2 . As was noted in ||16|| . 
deviations of the field D from the vacuum decrease with increasing mass m D , as well as 
with increasing number of colours N c , and the method works for any m D (and also for 
charm and bottom quantum numbers). 



The excitation frequencies u and w are: 



N C B 
86 



F,B 



F,B 



UJF,B 



N C B 
86 



F.B 



{HF,B + 1) 



(15) 



As was observed in flI7|, the difference uj f ,b -uf,b = N c B/(AQ f ,b) coincides, to the leading 



order in N c , with the expression obtained in the collective coordinate approach [|18 
At large m D the fi FjB — 4:fh D (r B <d FiB ) 1/2 /(N c B) and for the difference uj fa -u> f>b we o 
(N c = 3): 
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F,l 



1/2 



6 



F.B 
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B 



6 



F.B 



e 



FA 



3tain 



(16) 



Obviously, at large m D , the first term in (16) dominates and is positive if Ti/6f,i > T B /Q F ^ B . 
This is confirmed by looking at Table 1. Note also that the bracket in the first term in 
(16) does not depend on the parameters of the model if the background SU{2) soliton 
is calculated in the chirally symmetrical limit: both r and 6 scale like ~ l/(F 7r e 3 ). In a 
realistic case when the physical pion mass is included in (3) there is some weak dependence 
on the parameters of the model. 

The FSB in the flavour decay constants, i.e. the fact that Fk/F^ ~ 1.22 and F D jF^ = 
1.7±0.2, should be taken into account. In the Skyrme model this fact leads to the increase 
of the flavour excitation frequencies which changes the spectra of flavoured (c, b) baryons 
and puts them in a better agreement with the data [^0], |21[] . It also leads to some changes 
of the total binding energies of BS [I7|. This is partly due to the large contribution of 
the Skyrme term to the flavour moment of inertia Q F . Note, that in []16|] the FSB in 
strangeness decay constant was not taken into account, and this led to underestimation 
of the strangeness excitation energies. Heavy flavours (c, b) have not been considered in 
these papers. 

The terms of the order of N~ x in the Hamiltonian, which depend on the angular 
velocities of rotations in the isospin and the usual space and which describe the zero-mode 
contributions, are not crucial but important for the numerical estimates of the spectra 
of baryonic systems. To calculate them one should first obtain the Lagrangian of BS 
including all the terms upto 0(1/ N c ). The Lagrangian can be written in a compact form 
as: 

L ~ -M cl + 46 F , S [i>t£>(l - j) - \ipib&D - {rtDf - (D^D) 2 )] + 2Q F , B (W) + - P?~ 

-[T B ml + (f* - Fl)f B ]D*D(1 + - t) { tf D DW) ^fu)D^D, (17a) 

where d 2 = 2D^D and 

p = -i(b^fD - D^fb). 

For the axially symmetrical configurations, like the B = 2 torus, the term 6 TiB (w 3 - /3 3 ) 2 /2 
in (17a) should be substituted by 



sl = ^M( W3 - n n 3 - p 3 ) 2 + ^p(n? + nl), 



(176) 



where Qj are the components of the angular velocities of rotation in the usual space, 
OinOkn = e*;0/- Taking into account the terms — 1/N C we find that the canonical variable 
n conjugate to D is 



U = = ±® FjB [b{l~)-\rtbD+ itfD D]+i(Q T , B -20 F , B )uTD-iG TtB pTD+i^(l~ 



(18) 



^From (17a) the body-fixed isospin operator is: 



P f = dL/du) = 9 t ,bw + (26 F , S - @ t< b)P 



N r B 



D^tD. 



Using the identities 
and 



0tD = 2D^DD - (D r D + DW)D 



(3 2 = AD^DD^D - (D^D + D ] D) 



(19) 



(20a) 



(206) 

we find that the ~ 1/N C zero mode quantum corrections to the energies of skyrmions can 
be estimated [|15|, |16H as: 



AE 1/Nc 



26 



T,B 



■ [c F<B T r (T r + 1) + (1 - c F , B )I(I + 1) + {cf.b - c f ,b)If{If + 1) 



(21a) 



where / = I bf is the value of the isospin of the baryon or BS, which can be written also 
as: 



pf = 6 



1- 



28^ 



N c BO T 
46, 



dWd 



(22) 



with the operator I F = (tfrb - a T fa tT )/2. 

T r is the quantity analogous to the "right" isospin T r , in the collective coordinate 
approach [18]], and T r = pf — f F . The hyperfine structure constants c F . B and cf,s are 
defined by relations: 



1 - cf.b 
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(23) 



To take into account the usual space rotations the J-dependent terms should be added 
to (21a). For the axially-symmetrical configurations, like the B = 2 torus, they are equal 
to 0, [[J: 



AE 1/Nc = 



1 



2n 2 B 3 .B 



(J: 



6A2 , 



26 



J,B 



(21b) 



with Qj : b being the moment of inertia corresponding to the usual space rotations - or- 
bital moment of inertia, which is known to increase with increasing i?-number almost 

For such configurations the body-fixed 3-d component of 



proportionally to B 2 [|T7| , 

the angular momentum J^* and the nonstrange part of the 3-d component of the isospin 



bf 



111' 



bf 



r.3 



see, e.g. 



18, M and ref- 



(also body-fixed) are connected by the relation X 
erences therein). Realistic cases of multiskyrmions are intermediate between the case of 
incoherence of usual space and isospace rotations and the complete coherence, as in (216) 
for the rotation relative to the axis of axial symmetry. However, the J-dependent terms 
of the type (216) are cancelled mostly in the differences of energies of states which belong 
to the same SU(3) multiplet, i.e. when they have the same values of J, (p, q) and T ri3 . 

In the case of antiflavour excitations we obtain the same formulas, with the sub- 
stitution uj — * uj and — > -\x in (23). For example, 



6 



T,B 



-■F.B 



6 



F.B^F.B 



-(/< 



F,B 



(24) 



The excitation energies for antiflavours are close to — 0.59Gev for antistrangeness, 
— 1.75 Gev for anticharm and to ~ 4.95 Gev for antibottom. However, these numbers should 
be considered as lower bounds only since to calculate them we have used a simplified 
version of the bound state soliton model. 



4 Estimates of the spectra of multibaryons with strangeness, 
charm or bottom 



In the bound state soliton model, and in its rigid oscillator version as well, the states 
predicted do not correspond originally to the definite SU(3) or Sf7(4) representations. How 
this can be remedied was shown in []16[ ; see also Eq. (26) - (29) below. The quantization 
condition (p + 2q)/3 = B [[[(J, for arbitrary N c , changes to (p + 2q) = N C B + 3n qq -, where n qq 



is the number of additional valent quark-antiquark pairs present in the quantized states 



Tq| . For example, the state with B = 1, \F\ = 1, I = and n qq = should belong to the 
octet of (u,d,s), or (u,d,c), SU(3) group, if N c = 3; see also |jIB| . The state with B = 2, \F\ = 2 
and 7 = should belong to the 27-plet of the corresponding group, etc. The states having 
antiflavour quantum number, i.e. positive strangeness or bottom quantum number or 
negative charm should have the number of additional quark-antiquark pairs n qq > \F\ WB 



If 6f — > oo, Eqs. (21) go over into the expression obtained within the collective coordinate 
approach |10|, [17]]. In a realistic case, with Q T /@y ~ 2 - 2.7, the structure of (21) is more 



complicated. 

First we consider quantized states of BS which belong to the lowest possible SU(3) 
irreps (p,q) for each value of the baryon number, p+2q = 3B: p = 0, q = 3B/2 for even B, and 
p = 1, q= (3B- 1)/2 for odd B. For £ = 3, 5 and 7 they are 35, 80 and 143-plets, for B = 2, 4, 6 
and 8 - 10, 2~8, 55 and 9~l-plets. Since we are interested in the lowest energy states, we 
discuss here the baryonic systems with the lowest allowed angular momentum, ie J = 0, 
for B = 2, 4, 6 and 8. For odd B the quantization of BS encounters some difficulties (see 
p3|j), but the correction to the energy of quantized states due to the nonzero angular 
momentum is small and decreases with increasing B since the corresponding moment of 



inertia increases proportionally to — B 2 [22, 23 . Moreover, the J-dependent correction to 



the energy cancels in the differences of energies of flavoured and flavourless states. 

For the energy difference between the state with flavour F belonging to the (p,q) 
irrep, and the ground state with F = and the same angular momentum and (p, q) we 
obtain: 

AE B . F = |*V Fffl + [1(1 + 1) - T A T r + 1)] + (^-1)(^,b-2) Jf(Jf + 

where T r = p/2 and usually If = I - T r . Note that the moment of inertia 9t does not 
enter the difference of energies (25). Obviously, for "minimal" BS, i.e. those which do not 
contain additional quark-antiquark pairs, 

T r < 3B/2. (26) 

The maximal isospin carried by |F| flavoured mesons bound by (u,d) solitons satisfies 
another obvious relation: 

If = \F\/2. (27) 

Simple arguments allow us also to get the following restrictions on the total isospin of 
BS: 

\T r - \F\/2\ <I<T r + \F\/2 (28) 

and 

I < (3B- \F\)/2. (29) 

The lowest of the two upper bounds should be taken as the final upper bound. It is easy 
to check that our bounds correspond to the known SU(3) multiplets for each value of T r . 



For the B = 1 case, the difference of masses within the octet of baryons, Af and 
nucleon, T, F and A F , is 

\i\,t 3(1 -cf,i) 3(/x f ,i - 1) (1 - c F ,i) M-F.i - 1 , Qn N 

AAZApjV = — =WF,1 9 , AMs p A p = ~ = 7; ~ • (30) 



89 T ,i 



8(j? f1 Qf,i ' 



e T ,i 



2/UF,l©F,J 



Clearly, there are cancellations in Eq. (25) - the binding energy differencies of 
multiskyrmions. For states with maximal isospin I = T r + \F\/2 the energy difference can 
be simplified to: 



AE 



B,F 



\F\ 



uf,b + T r 



Mf,b - 1 



|F|+2) (^,b-1) 2 

y-F.B 



4(J>F>B&F,B 80F.B 

For even B T r = 0, for odd J3 we should take T r = 1/2 for the lowest ST/ (3) irreps. 



(31) 





Ae s= _i 


Ae c= i 


Ae 6= _! 


Ae s= _ 2 


Ae c=2 


Ae fc= _ 2 


Ae s= _ 3 


Ae c=3 


Ae s= _ 4 


2 


-0.047 


-0.027 


0.02 


-0.115 


-0.088 


0.02 


-0.205 


-0.183 


-0.316 


3 


-0.042 


-0.010 


0.04 


-0.098 


-0.040 


0.06 


-0.168 


-0.064 


-0.252 


4 


-0.020 


0.019 


0.06 


-0.051 


0.022 


0.10 


-0.092 


0.013 


-0.144 


5 


-0.027 


0.006 


0.05 


-0.063 


0.001 


0.08 


-0.108 


0.019 


-0.160 


6 


-0.019 


0.016 


0.05 


-0.045 


0.023 


0.10 


-0.078 


0.028 


-0.117 


7 


-0.016 


0.021 


0.06 


-0.041 


0.033 


0.11 


-0.070 


0.037 


-0.105 


8 


-0.017 


0.014 


0.02 


-0.040 


0.021 


0.03 


-0.068 


0.020 


-0.100 



Table 2. The binding energy differences Ae s , c ,b are the changes of binding energies of lowest BS with 
flavour s, c or b and isospin I — T r + |-F|/2 in comparison with usual u, d nuclei, for the flavour numbers 
S = -1, -2, -3 and -4, c = 1, 2 and 3, b = -1 and -2 (see Eq. (32)). 



It follows from (30) and (31) that when some nucleons are replaced by flavoured 
hyperons in BS the binding energy of the system changes by 



Ae_e p = \F\ 



3(mf,i ~ 1) T VF.B-l (1*1+2) (mf.b-1) 2 



8fi F1 &F,i 4/xf,b6fb 86 f ,b V f ,b 



(32) 



For strangeness Eq. (32) is negative indicating that stranglets should have binding energies 
smaller than those of nuclei, or can be unbound. Since Qf.b and 6 t ,b increase with 
increasing B and mo this leads to the increase of binding with increasing B and mass of 
the flavoured state, in agreement with ||17|| . For charm and bottom Eq. (32) is positive 
for B > 3 or 4. It follows from Table 2 that dibaryons with strangeness or charm quantum 
number are probably unbound, but those with b = -1 or b = -2 could be bound. The 
multibaryons with B > 4 and S = -1 can be bound, as well as multibaryons with c = 1, 2 
or 3, or bottom b = -1, -2. 

Had the momenta of inertia of BS at small values of B been proportional to the 
baryon number B, then the values of fi, excitation frequencies uj f and coefficients c would 
not have depended on B at all. In this case the binding energy would have consisted only 
of its classical part and a contribution from zero modes; the difference of w's would have 
been absent in this case. 

The nuclear fragments with sufficiently large values of strangeness (or bottom) 
may have been found in experiments as fragments with negative charge Q, according to 
the well known relation, Q = T 3 + (B+S)/2 (similarly for the bottom number). Recently one 
event of a long lived nuclear fragment with mass about lAGev was reported in [[MJ . Using 
the above formulas it is not difficult to establish that this fragment may be the state 
with B = —S = 6, or B = 7 and strangeness S = -4, or -3, see also Table 3 below. Greater 



values of strangeness are not excluded since the method used here overestimates the 
flavour excitation energies, especially for smaller baryon numbers and for the strangeness 
quantum number. 

Another case of interest involves considering the BS with isospin 1 = 0. In this case 
I F = T r = \F\/2 and so such states do not belong to the lowest possible SU(3) multiplet for 
each value of B (except for the case \F\ = l). For the energy difference between this state 
and a flavourless state belonging to the same SU (3) irrep it is easy to obtain: 



AE 



B.F 



\F\ 



\F\+2) {hf,b-1) 



86 



F,B 



(33) 



For the difference of binding energies of such a state and the ground (u, d) state with 
lowest values (p" 



, m ,<j mm ) we have the following estimate: 



Acb,f = \F\ 



where TL 



U>FA ~ I^F,B ' 



3(/x Fjl -l) , (|F|+2)( MFjS -l) 



86 



FA^p.l 



86 



F.B 



^F,B 



26 



T.B 



-[|F|(|F| +2)/4 -T™ n (T; 



HI)], 
(34) 

0, or 1/2. Using this formula we find the values given in Table 3. For 



example, the B = 2, |F| = 2 state discussed previously in |T8[ and later in [|16| belongs to 
the 27-plet of the corresponding SU (3) group. In the case of strangeness it has appeared 
already, probably, as a virtual level in the AA system pi . 
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Ae s= _i 


Ae c= i 


Ae 6= _i 


Ae s= _ 2 


Ae c=2 


Ae b= _ 2 


Ae s= _ 3 


Ae c=3 


Ae 6= _ 3 


2 








-0.075 


-0.029 


0.02 








3 


0.000 


0.034 


0.07 








-0.083 


0.002 


0.09 


4 








-0.047 


0.030 


0.09 








5 


-0.003 


0.032 


0.06 








-0.060 


0.035 


0.12 


6 








-0.044 


0.025 


0.09 








7 


0.000 


0.040 


0.07 








-0.042 


0.068 


0.15 


8 








-0.039 


0.023 


0.03 









Table 3. The binding energies differences of lowest flavoured BS with isospin 1 = and the ground state 
with the same value of B and I = or / = 1/2, see Eq. (34). The first 3 columns are for |F| = 1, the next 
3 columns - for \F\ = 2, and the last 3 - for |F| = 3. The state with the value of flavour |F| belongs to the 
SU(3) multiplet with T r = |F|/2. 



We can see from this Table that the B = 7, S = -3 state has a binding energy smaller 
than the (u, d) nucleus by 42 MeV, i.e. it can still be stable with respect to the strong decay, 
if we take into account the uncertainty of our estimates (recall that the nucleus 7 Li has 
the total binding energy 39 Mew.) The state with isospin equal to 2 - maximal value for 
S = -3 within the (1, 10) 5(7(3) multiplet - has somewhat greater binding energy, see Table 
2 . The difference of energies of states with isospin / = I max = T r + \F\/2 and 1 = and the 
same value of F can be written as 

e- - e ,=0 = - ACT = HM±1> 

8 \<C>F,B[J.F, 

At large \F\ this is approximately given by 

^-^=^-^-(l-2c F , fl ) (356) 

At large B and F the isoscalar states have smaller energy if cf,b < 0.5, see also Table 1. 

It is of interest to consider the case corresponding to the bulk of "flavoured matter" , 
i.e. p = q = B = \F\. Such "multilambda" states with isospin equal to zero have the 



_±_\_T r (T r + l) +TAFl ,F, B ~1 

6t,b / 26t,b 46f,s Hf,b 

(35a) 



following value of Ae for B > 1: 



Ae~ \F\ 



F\ + 2 / jJ-F.B - 1 



3(^f,i - 1) 



(36) 



At large \F\ the sign of this expression depends on the sign of the difference (hf.b - 
l)/(©F,BA < lB) _i /®T,s- To draw the final conclusion we need the knowledge of the behaviour 
of the ratio 9 t ,b/6f,b at large B. For heavy flavours, c and b, ^f,b > 1 and second term 
in (35) is negative, unless 0t,b ~ hf,b<3>f,b which is not realistic (we have usually n s — 3, 
fi c — 15 and fi b ~ 73). So, for heavy flavours it is not possible to obtain, in this way, 
the bulk of flavoured matter as quantized coherent multiskyrmions. Other possibilities 
remain to be investigated, e.g. flavoured skyrmion crystals. 

As in the B = 1 case [26], the absolute values of masses of multiskyrmions are 
controlled by the poorly known loop corrections to the classic mass, or the Casimir 

the renormalization procedure is 



18 



energy. And as has been done for the B = 2 states, 
necessary to obtain physically reasonable values of these masses. As the binding energy 
of the deuteron is 30 MeV instead of the measured value 2.2 MeV we see that ~ 30 MeV 



characterises the uncertainty of our approach JlT], [18]]. But this uncertainty cancels in 



the differences of binding energies calculated in Tables 2,3. 



5 Conclusions 

Using rational map ansaetze as starting configurations we have calculated the static char- 
acteristics of bound skyrmions with baryon numbers up to 8. The excitation frequencies 
for different flavours - strangeness, charm and bottom - have been estimated using a rigid 
oscillator version of the bound state approach of the chiral soliton models. One notes 
that, in comparison with strangeness, this approach works even better for c and b flavours 



BQ, 21]. Our previous conclusion that BS with charm and bottom have more chance to 



be bound by strong interactions than strange BS [|17| is reinforced by the present inves- 
tigation. Estimates of the binding energy differences of flavoured and flavourless states 
have some uncertainty, about few tens of Mev, but the tendency for charm and bottom 
to be bound stronger than strangeness is very clear. 

A natural question now arises as to how these results depend on the choice of the 
parameters of the model. A set of parameters, which has been used extensively in the 
literature, is, e.g. the set introduced in || where the masses of the nucleon and the A 
isobar have been fitted in the massless case and with physical pion mass, correspondingly. 
In view of the large negative contribution of the loop corrections, or the Casimir energy, 
we feel that this choice of parameters cannot be taken too seriously. But calculations 
show that our results hold for this choice too. The energies of the flavour excitation 
are somewhat smaller, however: for example, for strangeness uj s = 255 MeV for 5 = 1 and 
249 MeV for B = 4, and we take F K /F n = 1.22 {F n = 108 Mew, e = 4.84). Similar changes take 
place for charm and bottom, and the conclusion that charmed or bottomed BS have good 



chances to be stable against strong interactions remains valid |fl7 |. 



It should be kept in mind that corrections of the order of 1/N% can lead to some 



change of our results. For example, the flavour moment of inertia changes to [|T8 



6 F - 6^ 0) - D*D F ° g F " J(l- C/ )(2 - c f )d 3 r, (37) 

Decrease of the moment of inertia could lead to some increase of the zero-mode quantum 
corrections. 



The apparent drawback of our approach is that the motion of the system into 
the "strange", "charm" or "bottom" directions is considered independently from other 
motions. Consideration of the BS with "mixed" flavours is possible in principle, but its 
treatment would be more involved (see, e.g. [27] where the collective coordinate approach 
to the quantization of SU(n) skyrmions has been investigated). 

Our results agree qualitatively with the results of |28| where the strangeness exci- 
tation frequencies had been calculated within the bound state approach. The difference 
is, however, in the behaviour of excitation frequencies: we have found that they decrease 
when the baryon number increases from B = 1 to 4, thus increasing the binding energy of 
the corresponding BS. This behaviour seems to be quite natural: there is an attraction 
between K, D or B meson field and B = 1 nucleon, and the attraction of a meson by 2, 3 
etc. nucleons is greater. Similar results hold for ordinary nuclei: the binding energy of 
a deuteron is 2.22 Mew only, for B = 3 it is about 8Mev, for B = 4 it is already 28 Mew, and 
saturation takes place soon. 

There is a further difference between the rigid oscillator variant of the bound 
state model we have used here and the collective coordinate approach of soliton models 
studied previously |10]-[]12|. In the collective coordinate approach involving zero modes 



of solitons with a rigid or a soft rotator variant of the model, the masses of baryons are 
usually considerably greater than in the bound state approach when the Casimir energies 



are not taken into account |26| , [29|| . One of the sources of this difference is the presence of a 
term of order N C /Q F in the zero-mode contribution to the rotation energy, which is absent 
in the bound state model. Recently it has been shown by Walliser, for the B = 1 sector 



within the SU(3) symmetrical (m K = m n ) variant of the Skyrme model [p9 l, that this large 



contribution is cancelled almost completely by the kaonic 1-loop correction to the zero- 
point Casimir energy which is of the same order of magnitude, N® [[29 . This correction 



has also been recently calculated within the bound state approach to the Skyrme model 
13C 



The charmed baryonic systems with B = 3, 4 were considered in |3T] within a 
potential approach. The B = 3 systems were found to be very near the threshold and the 
B = 4 system was found to be stable with respect to the strong decay, with a binding 
energy of ~ lOMeV. 

Experimental searches for the baryonic systems with flavour different from u and 
d could shed more light on the dynamics of heavy flavours in few-baryon systems. The 
negative charge fragment seen in the NA52 CERN experiment []25| could be explained in 
our approach as a quantized B = 7 skyrmion with strangeness S = -3 or -4. The other 
possibility is B = 6 and S = -6 or -7. The value of strangeness can be greater since 
the rigid oscillator version of the model we consider here overestimates the strangeness 
excitation energies. 

The threshold for the charm production on a free nucleon is about l2Gey, and 
for the double charm it is ~ 25.2 GeV. For bottom, the threshold on a nucleon is — 70 
GeV. However, for nuclei as targets the thresholds are much lower due to the two-step 
processes with mesons in intermediate states and due to the normal Fermi-motion of 
nucleons inside the target nucleus (see, e.g. |32j). Therefore, the production of baryons 
or baryonic systems with charm and bottom should be feasable in proton accelerators 
with energies of several tens of GeV, as well as in heavy ions collisions. 

Let us finish by adding that a shortened and less complete version of this paper 
is available ||33|| . The results obtained recently in [33] within the detailed version of the 



bound state approach are in fair agreement with our, but the binding obtained in [^4 



is smaller than what we have found here. It should be noted that, in difference from 
3^|, we have used the empirical values of flavour decay constants, taken into account the 



l/N c zero modes contributions to the energy of multibaryons and have considered only 



the difference of binding energies of flavoured BS and the ground states where many of 
uncertainties cancel out. 

We are indebted to B.E. Stern for help in numerical calculations and to J.Madsen 



for informing us about the result of NA52 experiment ||25| . This work has been supported 

by the UK PPARC grant: PPA/V/S/1999/00004. 
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